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ABSTRACT 

Context. The search for extra-solar planets similar to Earth is becoming a reality, but as the level of the measured radial-velocity 
reaches the sub-ms~\ stellar intrinsic sources of noise capable of hiding the signal of these planets from scrutiny become more 
important. 

Aims. Other stars are known to have magnetic cycles similar to that of the Sun. The relationship between these activity variations and 
the observed radial-velocity is still not satisfactorily understood. Following our previous work, which studied the correlation between 
activity cycles and long-term velocity variations for K dwarfs, we now expand it to the lower end of the main sequence. In this first 
paper our aim is to assess the long-term activity variations in the low end of the main sequence, having in mind a planetary search 
perspective. 

Methods. We used a sample of 30 M0-M5.5 stars from the HARPS M-dwarf planet search program with a median timespan of 
observations of 5.2 years. We computed chromospheric activity indicators based on the Can H and K, Ha, Hei D3, and Nai Dl and 
D2 lines. All data were binned in to average out undesired eff'ects such as rotationally modulated atmospheric inhomogeneities. We 
searched for long-term variability of each index and determined the correlations between them. 

Results. While the ^Scaii, Ha, and Nai indices showed significant variability for a fraction of our stellar sample (39%, 33%, and 
37%, respectively), only 10% of our stars presented significant variability in the Hei index. We therefore conclude that this index is 
a poor activity indicator at least for this type of stars. Although the Ha shows good correlation with *Scaii for the most active stars, 
the correlation is lost when the activity level decreases. This result appears to indicate that the Caii-Hof correlation is dependent on 
the activity level of the star. The Nai lines correlate very well with the Scaii index for the stars with low activity levels we used, and 
are thus a good chromospheric activity proxy for early-M dwarfs. We therefore strongly recommend the use of the Na i activity index 
because the signal-to-noise ratio in the sodium lines spectral region is always higher than for the calcium lines. 

Key words. Techniques: spectroscopic - Stars: activity - Stars: late-type 

1, Introduction duced by extra-solar plane ts (iQueloz et al.1 120011: iBonfils et aP 

'< ^ . 120071: iHuelamo et al.ll2008l) . It is therefore extremely important 

. ^ The mcrease m precision of the radial-velocity mstruments is ^^^^^ understand how stellar activity produce these 

k> leading to the detection of smaller reflex semi-amplitude signals ^^^^^^^ ^^^^ corrected for. 

induced by extra- solar plan ets (e.g. the 1.9 M© planet discov- 
er ; ered bv iMavor etani2009h . Also, as the observation timespan Studies conducted at the Mount Wilson observatory uncov- 
increases, the detection of planets in longer orbital periods is be- ered that many solar-like stars have magnetic cycles similar to 
coming possible (e.g. Wrig ht et al. 2008). that of the Sun (Wilson 1978; Baliunas et al. 1995). These stud- 
As these keplerian semi-amplitudes reach the sub-ms'^ ies concluded that around 85% of the stars have stellar activity 
level, the signals induced by intrinsic sources (e.g. oscillations, variability (|Baliunas et al.||l998|). Although the Mt. Wilson sur- 
rotating spots and plages, inhibition of convection) become vey covered a broad range of spectral types, only one M-dwarf 
more significant (e.g . |Saar & D onahue 1997: Sa ntos etal.l|2Q^ was included in the sample. This star, HD 95735, has an aver- 
Paulson et al.l l200a iBouc^^ 2004: Meunier et al.l |20TqI: age S mw value of 0.424 and a variable activity cycle (with the 
Dumusque et alTl2011aabl: iBoisse et al.ll2011h . It is weU-known cycle period poorly defined as yet). More rece ntly, evidence for 



that steUar activity can be a source of radial-velocity "noise" cyclic activity was found for a few M dwarfs. ICmcunegui et al.| 

and can even induce periodic RV signals similar to those pro- UQQM claimed a magnetic cycle with a -442 day period for 

the dMe 5.5 Prox Centauri and iDiaz et al.1 (l2007bl) found a P 

^ Based on observations made with the HARPS instrument on the = ~763 day activity period for the dMe 3.5 spectroscopic bi- 

ESO 3.6-m telescope at La Silla Observatory under programme ID nary G1375. Aditionally, activity cycles were also claimed to be 

072.C-0488(E) present on G1229A (Ml/2) and G1752A (M2.5) with periods of 
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~4 yr and ~7 yr, respectively (iBuccino et al.ll201Ql) . Apart from 
these few cases, not much is known about the magnetic cycles 
of the largest population of stars in the solar neighborhood. 

In a recent paper we studied a sample of seven late-G and 
early-K dwarfs with well-known magnetic cycles and compared 
the activity level of these stars w ith their radial-velo city varia- 
tions over a timespan of five years (ISantos et al.l201Ql) . We found 
that generally the S index was correlated with Ha and anti- 
correlated with the He i index. Also, the long-term variations of 
the S index could be detected in the cross-correlation function 
parameter's line bisector span, full-width-at-half-maximum and 
contrast, implying that these parameters could be used to follow 
the activity cycles of these stars. 

In this first paper we extend this study to the lower end of the 
main sequence and use a more extended sample. We study for 
the first time the eff'ect of long-term chromospheric activity of a 
sample of dwarfs of spectral types from MO to M5.5, covering 
the range of stars from partially convective to the beginning of 
fully convective interiors. We used the same HARPS M-dwarf 
sample as in the search of extra- solar planets in the southern 
hemisphere (see Bonfils et al. 201 1, submitted). 

The outline of this paper is as follows: we first describe our 
sample and observations in Section^ in Section[3]we explain the 
computation of our four activity indices, we search for activity 
variability in Section Hand for maxima or minima of activity in 
Section[5] compare the indices in Section^ and our conclusions 
are finally presented in Section [71 



2. Sample and observations 

The sample comes from the HARPS (iMavor et al.l l2003h M- 
dwarf planet search program, which corresponds to a vol- 
ume limited selection of stars brighter than V = 14 mag and 
with a projected ro tational velocity vsin/ < 6.5 kms"^ (see 
iBonfils et al.l l201l[ submitted). They obtained high-resolution 
spectra (spectral resolution = 115 000) for these stars from 2003 
to 2009, with 15 min integration times. 

For the calculation of our S inde^H we selected all spectra 
with a signal-to-noise ratio at spectral order 6 (corresponding to 
the Can K line, -3933 A) higher than 2, the value at which the 
relation between the index and S/N disappears. All other indices 
were not aff'ected by this selection except for cases where we 
compared them with 5caii- The data were then nightly averaged 
and bins of 150 days were used to suppress any rotationally mod- 
ulated signals caused by short-term activity. Each bin included 
at least three nights, and cases where this condition was not met 
were discarded. At this stage only stars with four or more bins 
were considered for the rest of this study. This selection enabled 
us to compare the diff'erent parameters and look for correlations 
with enough data points. From now on, we refer the binned data 
unless specified. We ended up with a sample of 30 stars matching 
these conditions (23 for the S index) with spectral types ranging 
from MOV to M5.5V. 

The errors used for each bin are the errors on the average, 
cr/ ^fN, where cr is the root-mean- square (rms) of the nightly 
averaged data in each bin and N the number of measurements 
used to calculate the average. 

Table [T] presents some basic parameters and the observation 
log for these stars. The start and end dates, number of nights. 



^ From now on the index based on the Ca ii lines will be referenced 
as S index or »Scaii- 



timespan of observations, and signal-to-noise ratios were calcu- 
lated for the nightly averaged data without using the S /N > 2 
selection. The timespan of observations ranges from 1.9 to 6.5 
years, which should be sufficient to detect activity cycle varia- 
tions in these stars if they are present. 

3. Activity indices 

The Can H and K lines are well-known and widely used stel- 
lar mag netic activity prox ies. Our Can index was computed as 
in Boiss e et al.l (l2009h and lSantos eTaP (HOOO) with sHght mod- 
ifications to the wavelength bands at the core of the lines. We 
chose bands of 0.6 A centered at the Can H (/139 68.47) and K 
(/1393 3.66) lines instead of the 1 .09 A band used by lBoisse et al.l 
(l2009h because this last one includes parts of the line wings 
and thus undesirable photospheric flux (see Fi gure [H a similar 
smalle r band for the K line was also used by [Livingston et al.] 
(l2007l) for the same reasons). The sum of the flux in the two H 
and K lines was then weighted by the square of their respective 
errors taken as Va^ where N is the number of counts inside each 
band. Our reference bands were two 20 A windows centered at 
3900 and 4000 A. Because the Mt. Wilson survey only includes 
one M dwarf, we did not have enough stars to calibrate our 5'caii 
index to the S mw scale. Furt hermore, we could not compute a 
photospheric free 7?^^ index (iNoyes et al.lll984l) because there 
is no calibration for stars with B-V higher than 1.2. Acompar- 
ison of the indices with color is discussed in Section [6^ 

The Ha index is also a widely used chromospheric indicator, 
and some authors recommend it for the study of later-type stars 
owing to the in creasing flux in the red part of the spectra for 
these s tars (e.g. iPasquini & Pallavicinil [19911: ICincunegui et al 
l2007bh . Our Ha index was also computed as in iBoisse et al 
(l2009l) but using a slightly broader central band. Instead of the 
0.6 A window centered at 6562.808 A, we used a 1 .6 A band that 
seems to inclu de more contribution from chro mospheric activity 
(see FigurefTt. lPasquini & Pallavicinil (Il99ll) obtained bands in 
the region 1.6-1.7 A for the chromospheric contribution of the 
Ha line by subtracting the spectra of inactive stars to those of 
active st ars of the sam e spectral type (the same band was also 
found bv lHerbigll985l) . Recently. [Cincunegui et all (l2007bh also 
used a window of 1.5 A to calculate the flux in the Ha line. We 
divided the flux in the central line by the flux in two reference 
lines of 10.75 and 8.75 A centered at 6550.87 and 6580.31 A, 
res pectively. 

iDiaz et al.l (l2007bh proposed that the Nai Dl and D2 lines 
could also be used to follow the chromospheric activity level of 
very active late-type stars with (B-V) > 1.1. In this sample 
we do not have very active stars, but we decided to study these 
lines as a complement to the other widely used activity prox- 
ies. The Nai D lines provide inform ation about th e conditions in 
the middle-to-lower chromosphere (lMauasll2000l) and thus are a 
good complement to the upper chromosphere indicator Ha and 
lower chromosphere proxy Can index. We first com puted our 
Na I index in a similar fashion as iDiaz et al ] (I2007bl) by calcu- 
lating the average flux in the core of the Dl (/15895.92) and D2 
(/15889.95) lines using 1 A bands. We used two reference bands 
with windows of 10 and 20 A centered at 5805.0 and 6090.0 A, 
respectively. For each band we selected the 10 higher flux values 
and calculated the average. The average values of the two Dl and 
D2 lines were then divided by the average of the two reference 
bands. We noticed that the 1 A bands include a significant part of 
the line wings and thus possible photospheric contribution to the 
flux. Therefore, we tried narrower bands of 0.5 A for the lines 
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Table 1. Basic parameters of the stellar sample and observation log. 



Star 


Sp. Type^ 


V 
[mag] 


V-I^ 
[mag] 


BJD,,^,, - 2400000 
[days] 


BlDend - 2400000 
[days] 


^ nights 


T 

^ span 

[years] 


<SIN> 
(sp. order 6) 


GJ361 


M1.5V^ 


10.40 


2.01 


54455.84 


55169.84 


34 


2.0 


3.9 


GJ2049 


Ml 


11.17^ 


1.77 


54455.70 


55155.86 


25 


1.9 


3.2 


GJ3218 


M2Y' 


11.12 


1.65^ 


54396.89 


55141.67 


40 


2.0 


2.0 


Gil 


M3V 


8.57 


2.13 


52985.60 


55048.83 


43 


5.6 


11.5 


G1176 


M2.5V 


9.97 


2.25 


52986.71 


55132.82 


70 


5.9 


4.2 


G1205 


M1.5V 


7.92 


2.08 


52986.73 


54386.81 


75 


3.8 


15.0 


G1273 


M3.5V 


9.89 


2.71 


52986.77 


55126.87 


61 


5.9 


6.5 


G1382 


M2V 


9.26 


2.18 


52986.84 


54174.65 


30 


3.3 


8.0 


G1393 


M2V 


9.76^ 


2.26 


52986.86 


54883.74 


28 


5.2 


6.9 


G1433 


M2V 


9.79 


2.15 


52989.84 


55057.47 


60 


5.7 


5.3 


G1436 


M2.5 


10.68 


2.02 


53760.83 


54999.46 


105 


3.4 


2.6 


G1479 


M3V 


10.64 


1.90 


53158.55 


54571.70 


57 


3.9 


4.2 


G1526 


M1.5V 


8.46 


2.07 


53158.60 


55001.58 


34 


5.0 


11.6 


G1551 


M5.5V 


11.05 


3.62 


52684.86 


55057.54 


37 


6.5 


1.2 


G1581 


M2.5V 


10.57 


2.53 


53152.71 


55056.53 


128 


5.2 


3.6 


G1588 


M2.5V 


9.31 


2.40 


53152.75 


54956.80 


25 


4.9 


9.2 


G1667C 


M2V 


10.22 


2.08^ 


53158.76 


55053.69 


147 


5.2 


4.0 


G1674 


M3V 


9.36 


2.40 


53158.75 


54732.48 


44 


4.3 


8.4 


G1680 


M1.5V 


10.14 


2.27 


53159.71 


55057.65 


28 


5.2 


4.3 


G1699 


M4V 


9.54 


2.52 


54194.89 


55054.65 


25 


2.4 


6.6 


G1832 


Ml V 


8.67 


2.18 


52985.52 


55122.65 


57 


5.9 


10.4 


G1849 


M3V 


10.42 


2.50 


52990.54 


55122.62 


48 


5.8 


4.5 


G1876 


M3.5V 


10.17 


2.77 


52987.58 


55122.63 


37 


5.8 


4.8 


G1877 


M2.5 


10.55^ 


2.43 


52857.83 


54704.73 


39 


5.1 


2.8 


G1887 


M2V 


7.34 


2.02 


52985.57 


54392.68 


63 


3.9 


15.7 


G1908 


Ml V 


8.98 


2.04 


52986.58 


55057.85 


50 


5.7 


11.6 


HIP 12961 


MOV" 


9.7 


1.64 


52991.63 


55109.79 


45 


5.8 


3.4 


HIP 19394 


M3.5V^ 


11.81 


2.5 


52942.80 


55105.79 


29 


5.9 


1.2 


HIP38594 


M^ 


9.96 


1.64 


52989.79 


55105.90 


16 


5.8 


2.8 


H1P85647 


MOV^ 


9.59 


1.85 


53917.75 


55117.49 


27 


3.3 


3.8 



References. fa) lBonfils et"an ( l201lL submitted) unless individually specified; fbl lHawlev etaP (119961) : (c) lKoen et an(l201Ql) : (d) lESAl (119971) : (e) 
Simbad (http://simbad.u-strasbg.fr/simbad/) 

pie, G1479 with < ^caii > = 0.0960 but slightly diflTerent color 
(y - / = 1.90). The dashed lines mark the line centers, while the 
dotted lines delimitate the band window. The flux in all figures 
was normalized by dividing it by the mean flux in the reference 
bands of each index. 

For these four stars the Can and Nai lines show a simi- 
lar behavior. Both increase in emission with increasing activity 
level in the S -index. On the other hand, the Ha and He i lines 
do not show the same trend. The Hof first decreases its emis- 
sion as activity increases and then exceeds the level of the star 
with less activity, ultimately arriving at a level similar to that of 
the continuum (if the activity increases even more it is expected 
to become an emission line). This behavior was explained by 
[Cram&]V[ullanl(ll979h through non-LTE model chromospheres 
for dK and dlVL stars. These authors explained that the stronger 
absorption in the line with increasing chromosphere heating rate 
is caused by an increase of the n = 2 bound level population 
of hydrogen and that when the heating reaches a certain level, 
the chromosphere electron density becomes so high that the line 
becomes collisionally excited, which increases the emission at 
its core. The He i D3 triplet is observed in Figure [T] in absorp- 
tion. This is generally the case for surges, eruptive prominences, 
and weaker flares, whil e emission is correlated with more in- 
tense flares (IZirinlll988b . The behavior observed in the figure is 
that of an anti-correlation with Can: as the Can lines emission 
increases, there is a tendency for the He i lines to increase in ab- 
sorption. 



and noticed that the correlation with the 5'caii index improved 
drastically (see the results section). These narrower bands were 
formally adopted for the rest of this work . 

The He i D3 line (/15875. 62) was observed to show good spa- 
tial correlation with plages (lLandmanlll981) and can therefore 
be used as an activity proxy (see lSaar et all 19971: llMontes et al.l 
1 19971 and references therein). We used a 0.4 A band centered at 
the line and divided its flux by the flux in two reference windows 
of 5 A centered at 586 9.0 and 5881.0 A, following the procedure 
in lBoisse etan(l2009h . 

The errors of the four indices were estimated by diflTerenti- 
ating the respective equations, and taking into account the flux 
in each band. But because we are interested only on the long- 
term variations of the indices and we binned the data into 150 
day bins, the final errors are the statistical errors on the mean as 
explained in Section O The individual errors of the nightly ob- 
servations are always lower than the estimated statistical errors 
on the mean used for the binned data. 

Figure [T] shows spectra in the region around the lines used 
to calculate the activity proxies. The four lines in each plot 
represent three stars with diflTerent activity levels but similar 
color (y - / ~ 2.10) plus one of the most active stars in 
the sample. By increasing activity, the blue line stands for GU 
(< ^caii > = 0.0147, y - / = 2.13), the black line for G1526 
(< ^caii > = 0.0307, V - I = 2.07), the red line represents 
G1205 (< ^caii > = 0.0920, y - / = 2.08), and the green line 
illustrates the lines for one of the most active stars in the sam- 
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A [A] A [A] 



Fig. 1. Cores of the lines used to calculate the indices for three 
stars with similar color (y-/~2.10) and different activity levels 
plus one of the most active stars in the sample (green, G1479). In 
blue is GU, black is G1526, and red is G1205. 



4. Activity indices variability 

After obtaining the four activity indices for our sample we now 
try to infer which stars show long-term variability and if mag- 
netic activity cycles are present in the data. With this in mind we 
calculated the rms (o-^) and average errors (< cr/ >) for the in- 
dices and compared them using an F-statistics with the F-value 
computed a s F = crlj < p-j >^ following a similar procedure as 
in IZechmei ster et aP (l2009b . This approach will give the proba- 
bility, given an F-value, that the observed variability (cTe) can be 
explained by the mean errors of the data (ct/). In this case, the er- 
ror of the data depends on the high-frequency activity variations 
of the nightly averaged data and number of bins (as explained 
in Section O. Therefore, lower values of probability (P(F)) will 
indicate that the variability is not due to the scatter induced by 
the errors. The results of this study are presented in Table [2l For 
the S index we used the S/N > 2 selection, resulting in Nhins(S) 
bins. Where there are less than four bins, the columns relating 
to the variability of the 5'caii index appear empty. For the other 
indices we used no selection, resulting in Nt,ms shown in column 
3. The probability that the variation is caused by the errors is 
presented in the form of a P- value, where values lower than 0.05 
(95% significance) are highlighted in bold. 

Because we averaged the data by 150 day bins, we inten- 
tionally removed the short- timescale variations, and only long- 
term variations are to be detected. As a result of this, we expect 
to find variability in stars with longer periods of observation, 
where their activity cycle pattern becomes detectable. Indeed, of 
the four stars with a timespan shorter than three years only one 
(G1699) shows significant variability, and just in one of the four 
indices. Furthermore, from the 18 stars that show variability, 13 
have timespans longer than or equal to five years. Indeed, 72% 
of the stars with timespans equal to or longer than five years 
present significant variability in at least one activity index. This 



is expected because the ability to detect magnetic cycle activity 
variations should increase with timespan. 

Figures[2l[3l|4l[5l[6l and[7]show the time- series of the four ac- 
tivity indicators for the selected stars as having significant vari- 
ability. Small points without errorbars are the data averaged per 
night, points with error bars represent the binned data, where 
the errors are the standard error on the mean, the numbers are 
number of nights used in each bin, and the dashed lines are the 
second-order polynomials best-fitted to the data. In each plot 
there is a reference to the peak-to-peak variation (A) and to the 
standard deviations of the binned data (cr). The values in paren- 
thesis are the percentages of the variations relative to the average 
values. 

Although we may not have enough data to detect full cycles, 
there are some hints for the cycle phase at which some of these 
stars could be. For example, in Figure [3] the time-series of 5'caii 
for G1433 seems to show a maximum of activity for this star, 
while for G1581 a minimum is clearly visible. 
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Fig. 2. Variations of 5'caii, Ha , Na l and He l with time for the 12 stars with long-term variability (see Section (4]). For visualization 
purposes, the y-axis is defined as the mean of the data plus and minus 0.35 times the maximum of the average errors of each index. 
The X-axis is constant for all stars. Note that G1877 does not have data for the 5'caii index. 
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Fig. 3. Continued from Fig. [2l 
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Fig. 4. Continued from Fig. [3l 
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Fig. 5. Continued from Fig. |4l 
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Fig. 6. Continued from Fig. \5\ 
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Fig. 7. Continued from Fig. [6l 
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Table 2. Variability F-tests for the activity indices. Bold values indicate probabilities lower than 0.05 (95% significance level). The F-value used was F = crl/ < (Ti >^ where 
is the rms of the index (binned) and < cr/ > the average of the error on the mean for the binned activity data of each star. 
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4.1. 5'caii variability 

Although the spectral orders of the Ca ii H & K lines are the ones 
with the lower signal-to-noise ratio, the 5'caii index still shows 
a great number of stars with variability. From our sample of 23 
stars, 9 have P- values lower than or equal to 0.05, which rep- 
resents 39% of the sample. These stars are GU, G1273, G1433, 
G1581, G1588, G1667C, G1832, G1908, and HIP85647. If we in- 
clude the three stars with P- value below 0.1 (GJ361, G1876, and 
G1887; in general they have short observation timespans or only 
a few nights and thus any variations caused by their hypothetic 
cycle might be underestimated) we have 52% of the sample with 
significant mid-term activity variations. Moreover, all stars with 
significant variations (P < 0.05) in the 5'caii index also present 
significant variations in at least one of the other indices, which 
reinforces this variable behavior. 

4.2. Ua variability 

The Ha index indicates significant variability (P < 0.05) for 
10 stars (GU, G1433, G1436, G1581, G1588, G1699, G1849, 
G1877, HIP19394, and HIP38594) out of 30, representing 33% 
of the sample. All these stars, apart from G1699, HIP19394 and 
HIP38594, also have significant variations in the other indices. 
If we include the three stars with a P- values below 0.1 (GJ3218, 
G1273, and HIP85647), then around 43% of the sample shows 
variability in this index, and two of these stars present a strong 
variability in the ^caii and Nai indices. 



presence of magnetic cycles. We fitted straight lines and second- 
order polynomials to the time- series and calculated an F-test us- 
ing Fpoiy = (N- ^)(x\ope -x\oiy^lx\oir where;^^^^ is the chi- 
squared of the straight line model, X^^^iy the chi-squared of the 
polynomial model, and N the number of bins. The probability 
that the data are better fitted by a polynomial than by a straight 
line is given in Table [3] by low values of P{Fpoiy). We used P- 
values lower than 0.05 as a limit for statistically significant im- 
provements to the fit. 

The 5" Call index shows three stars having maxima or minima, 
namely GU, G1433, and G1581, whose data are statistically bet- 
ter fitted by a second-order polynomial than by a straight line. 
Both GU and G1433 show a maximum of activity in the time- 
series plots (Figures [2] and [3]), while G1581 presents a minimum 
(Figure 131). 

The Na l index shows the same behavior in the tests as the 
calcium lines apart from HIP85647. This star passed the test and 
shows a maximum in the time- series (Figure [7]) with a similar 
behavior also observed on 5'caii (P = 0.11) and Ha (P = 0.12). 
Two other stars passed the test for the Ha index, G1832 and 
HIP38594, but without similar variations observed for the other 
activity indicators. For G1581 we observe a maximum instead of 
the minimum produced by the Ca 1 1 and Na 1 1 lines. 

Two stars present significant maximum or minimum in He i. 
These are G1361, which has passed the test for the other three 
indices, and G1581 with a minimum. 



4.3. Nai variability 

In terms of variability, the Na i index seems to be one of the best 
we have. Out of 30 stars, 11 present significant scatter (G1273, 
G1433, G1436, G1526, G1588, G1667C, G1832, G1849, G1876, 
G1877, and G1908, representing 37% of the sample), and there 
are six more with 0.05 < P < 0.1 (GJ361, GJ2049, GJ3218, 
G1581, HIP38594, and HIP85647), which, if these are included, 
amounts to 57% of the sample. Apart from G1526, all stars with 
P < 0.05 also have a strong variability or hints of variability in 
the other indices. 



4.4. He I variability 

The He i index is the one with fewer stars presenting variability 
in the F-tests. Only three stars in the sample (G1205, G1908 and 
HIP85647, 10% of the sample) have P < 0.05 and five (G1176, 
G1551, G1680, G1887, and HIP12961) have 0.05 < P < 0.1, 
which yields an overall of 27% of the sample with variability. 
This indicates that this index, although in a region of the spec- 
trum with high signal-to-noise ratio, does not vary with a suffi- 
ciently high amplitude to be detected in M dwarfs. It is curious 
that G1205, although having a strong variability in the He i index, 
does not present any significant scatter either in 5'caii (P = 0-35) 
or Ha (P = 0.49). The same is true of G1176 and HIP12961. On 
the other hand, there are a large number of stars that have strong 
variations in the other indices, but their P- values for He i indicate 
no variability. Therefore, this index does not seem to be a good 
proxy of medium to long-term activity variations of early-M type 
dwarfs. 



5. Activity maxima and minima 

Apart from the variability tests we also performed F-tests to try 
to find maxima or minima in our data that could indicate the 



6. Comparision between activity indicators 

After assessing the variability for our four indices, we now com- 
pare them to assess their mid-term correlations. Table |4] shows 
the Pearson correlation coefficients between the activity indices 
for our sample (p). N^ins is the number of bins for each star 
(NbinsiS) for comparisons using the S index). The false-alarm 
probability (FAP) was computed by bootstraping the nightly av- 
eraged data, then binning the data every 150 days, and calculat- 
ing the coefficient for each of the 10 000 permutations. A signif- 
icant FAP was then chosen for values < 0.05 (95% significance 
level) and highlighted in bold. 



6.1. 5" Call vsUa 

Six stars (out of 23, corresponding to 26% of the sample) show a 
significant positive correlation between the 5'caii and Ha indices 
(Table 13]). Around 57% of the sample have a positive correlation 
between S and Ha higher than 0.5. Although it has been claimed 
that a strong correlation between these two indices exists, it is 
not clear if that is the case for all stellar types or levels of activity 



(e.g. 'Monies et al.lll995l:IStrassmeier et al.lll99Ql:lRobinson et al.1 



1990; Giamp apa et al 



iCincunegui et al.l 



1989h 



2007bh studied the correlation between 



these two indices for a sample of 109 stars ranging in spectral 
type from F6 to M5. These authors found a great variety of cor- 
relations between -1 and 1. 

More recently, IWalkowicz & Hawleyl (l2009l) compared si- 
multaneous spectra of Can and B aimer lines for a sample of 
M3 dwarfs and observed that the relationship between calcium 
and Ha is not linear: weak absorption of the Ha line can corre- 
spond either to weak or intermediate Ca ii K emission, as pro- 
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Table 3. F-tests for trends in our activity indices. The F- value used was Fpoiy = (N - '^Mxliope ~ ^j)oiy^^^j)oiy' "^^^^^ ^liope 
chi-squared of the straight line model, the chi-squared of the polynomial model, and N the number of bins. This test will give 
the probability that a second-order polynomial is a better fit to the data than a straight line. Bold values indicate probabilities lower 
than 0.05 (95% significance level). 



Star Nhins(S) Nmus Sq^h Ha Nai Hei 









P poly 


P(Fpoly) 


Fpoly 




Fpoly 


PiFpoly) 


P poly 


P(Fpoly) 


GJ361 


4 


4 


0.475 


0.72 


3.13 


0.37 


0.00785 


0.99 


414 


0.035 


GJ2049 


<4 


4 






0.187 


0.85 


0.0204 


0.98 


0.0733 


0.93 


GJ3218 


<4 


4 






5.44 


0.29 


5.46 


0.29 


0.00831 


0.99 


GU 


5 


5 


27.9 


0.019 


3.61 


0.12 


22.5 


0.023 


0.729 


0.34 


G1176 


6 


7 


1.33 


0.24 


0.579 


0.41 


0.442 


0.46 


0.827 


0.34 


G1205 


6 


6 


0.00533 


0.56 


0.0243 


0.56 


0.166 


0.53 


0.0851 


0.55 


G1273 


5 


6 


0.0998 


0.52 


0.0967 


0.55 


5.18 


0.059 


0.0761 


0.55 


G1382 


5 


5 


0.511 


0.39 


0.820 


0.32 


0.150 


0.50 


0.136 


0.51 


G1393 


<4 


4 






0.0149 


0.99 


27.6 


0.13 


0.0533 


0.95 


G1433 


5 


6 


17.6 


0.030 


5.15 


0.059 


18.4 


0.011 


0.0255 


0.56 


G1436 


6 


7 


0.196 


0.52 


2.79 


0.098 


0.769 


0.35 


0.111 


0.56 


G1479 


4 


4 


24.1 


0.14 


6.86 


0.26 


2.92 


0.38 


10.9 


0.21 


G1526 


4 


4 


6.02 


0.28 


0.666 


0.66 


9.94 


0.22 


0.168 


0.87 


G1551 


<4 


7 






0.193 


0.54 


0.795 


0.35 


0.510 


0.44 


G1581 


9 


10 


11.9 


0.0022 


12.2 


0.0010 


10.7 


0.0016 


4.45 


0.019 


G1588 


4 


4 


0.469 


0.72 


0.290 


0.80 


0.322 


0.78 


0.0387 


0.96 


G1667C 


8 


7 


0.713 


0.38 


3.55 


0.070 


1.26 


0.24 


4.06 


0.057 


G1674 


4 


5 


0.128 


0.89 


0.826 


0.32 


0.467 


0.40 


0.0445 


0.54 


G1680 


4 


4 


0.0706 


0.94 


1.48 


0.50 


4.78 


0.31 


0.905 


0.60 


G1699 


4 


4 


2.92 


0.38 


62.9 


0.089 


1.74 


0.47 


0.0186 


0.98 


G1832 


6 


5 


0.153 


0.54 


63.7 


0.0084 


0.0227 


0.54 


0.996 


0.29 


G1849 


6 


6 


0.285 


0.49 


0.673 


0.37 


0.325 


0.48 


1.08 


0.28 


G1876 


4 


4 


1.08 


0.56 


4.78 


0.31 


1.63 


0.48 


6.21 


0.27 


G1877 


<4 


5 






3.56 


0.12 


1.41 


0.24 


0.00635 


0.54 


G1887 


4 


4 


9.44 


0.22 


0.465 


0.72 


0.756 


0.63 


0.00450 


0.99 


G1908 


6 


8 


1.16 


0.27 


0.689 


0.39 


0.970 


0.30 


2.29 


0.11 


HIP12961 


6 


6 


0.0851 


0.55 


0.136 


0.54 


2.15 


0.16 


0.451 


0.43 


HIP19394 


<4 


6 






1.37 


0.23 


2.15 


0.16 


5.83 


0.051 


HIP38594 


<4 


4 






1399 


0.019 


37.3 


0.12 


0.322 


0.78 


HIP85647 


5 


5 


4.15 


0.11 


3.66 


0.12 


15.3 


0.034 


0.133 


0.51 



posed bylCram & Giampapal (Il987l) a nd previously observed b y 
IStauff'er & Hartmannl(ll986h (see also lRauscher & Marcvl2006h . 

Im eunier & Delfossd (l2009l) also studied these indices for the 
Sun, and although we used 150 day bins and they used individ- 
ual measurements, their discussion about the relation between 
Ca II and Ha can provide interesting insights here. These authors 
pointed out that if the timescale of the observations covered less 
than the solar cycle, the correlation could decrease to negative 
values. We do not have information about the length of the cy- 
cles of these stars (or if they have periodic cycles), but because 
our maxim um timespans are in the 6-ye ar range or less, the con- 
clusion by iMeunier & Delfossd (l2009l) for the low positive or 
negative values of the correlation coefficient might be a possi- 
bility. Another possibility to explain the low and negative values 
of the correlation coefficients might be explained by the differ- 
ent sensitivity of the indices to different activity phenomena. For 
the Sun, the surface coverage of plages as measured by the two 
indices is not exactly the same: it is smaller for Ha than for the 
Ca II lines (Meunier & Delfosse 2009). This will reduce the cor- 
relation between them. Furthermore, the H a line is more sensi- 
tive to the presence of dark filaments. As iMeunier & Delfossd 
(l2009l) pointed out, this will reduce the correlation coeflScient 
to values closer to zero if the filaments and plages are not shar- 
ing the same positions on the disk of the star. If the filaments 
are spatially well-correlated with plages, the correlation between 
the two indices will tend to negative values. Thus, the anti- 



correlation found for G1526 with a Pearson correlation coeflS- 
cient of -0.78 (FAP = 13%) can be an indication that the plages 
and filaments are spatially well-correlated for that star. Stars with 
strong positive correlations (and higher levels of activity) such as 
G1176, G1205, G1382, and G1479 probably have a saturated fil- 
ament contribution to Ha that therefore does not inffuence the 
index as much as the plage contribution. 

Figure [8] (left panel) shows the Pearson correlation coeffi- 
cient between these two indices against the mean values of the 
5" -index. Open squares are values with FAP < 0.01, open circles 
for 0.01 < FAP < 0.05, and "plus" ("-h") symbols for FAP values 
higher than 0.05. 

As was found by lCincunegui et all (l2007bl) , we obtain a great 
variety of correlations from around -0.8 to 1. Nevertheless, there 
is a tendency in our sample for the positive correlations. A 
trend can be observed when the correlation coefficient is plot- 
ted against the mean 5'caii^ for values of S lower than around 
0.035 there are no significant correlations and for S higher than 
-0.035, there are only positive correlations and some of them 
with statistical significance. This trend is not observed when the 
correlation is plotted against the mean values of Ha. 

An interpretation of this could be that because the Ha is 
more sensitive to filaments than the S -index, as the activity gets 
stronger (higher 5'caii) the contribution of plages becomes more 
important to the Ha index than the contribution coming from 
filaments, because their contribution saturates at a certain activ- 
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Table 4. Pearson correlation coefficients between the activity indices. FAPs calculated using bootstrap permutations (see text). In 
bold are given FAP values below 0.05. Nbins is the number of bins for each star. 



Star 






'5' Call 


vs Hof 


S Call 


vs Nai 


S Call 


vs He I 


Hof vs Na I 


Hof vs He I 


Na I vs He I 








P 


FAP 


P 


FAP 


P 


FAP 


p 


FAP 


p 


FAP 


p 


FAP 


GJ361 


4 


4 


0.58 


0.27 


0.73 


0.12 


0.43 


0.35 


-0.07 


0.46 


0.86 


0.13 


-0.24 


0.41 


GJ2049 


<4 


4 














-0.04 


0.47 


-0.89 


0.083 


0.43 


0.33 


GJ3218 


<4 


4 














1.00 


0.0006 


0.39 


0.36 


0.36 


0.37 


GU 


5 


5 


0.12 


0.41 


0.89 


0.031 


0.56 


0.24 


0.37 


0.33 


0.14 


0.45 


0.84 


0.067 


G1176 


6 


7 


0.85 


0.016 


0.33 


0.27 


0.76 


0.033 


0.40 


0.24 


0.77 


0.034 


-0.04 


0.48 


G1205 


6 


6 


0.98 


0.0015 


0.94 


0.0079 


0.58 


0.17 


0.83 


0.054 


0.57 


0.17 


0.69 


0.10 


G1273 


5 


6 


0.06 


0.46 


0.90 


0.024 


0.45 


0.25 


-0.40 


0.24 


0.80 


0.038 


0.08 


0.44 


G1382 


5 


5 


0.99 


0.0007 


0.98 


0.0017 


0.42 


0.24 


0.98 


0.0021 


0.46 


0.23 


0.56 


0.15 


G1393 


<4 


4 














-0.10 


0.45 


0.91 


0.044 


0.17 


0.42 


G1433 


5 


6 


-0.28 


0.32 


0.98 


0.0012 


-0.50 


0.29 


-0.32 


0.25 


-0.39 


0.30 


-0.31 


0.33 


G1436 


6 


7 


-0.06 


0.45 


0.87 


0.010 


0.42 


0.20 


0.49 


0.19 


-0.15 


0.38 


0.21 


0.32 


G1479 


4 


4 


0.99 


0.0044 


0.94 


0.034 


0.87 


0.060 


0.94 


0.033 


0.82 


0.10 


0.91 


0.050 


G1526 


4 


4 


-0.78 


0.13 


0.99 


0.013 


0.21 


0.36 


-0.86 


0.12 


0.42 


0.24 


0.10 


0.43 


G1551 


<4 


7 














0.59 


0.086 


0.96 


0.0007 


0.55 


0.11 


G1581 


9 


10 


-0.44 


0.094 


0.82 


0.0036 


0.75 


0.024 


-0.75 


0.0035 


-0.14 


0.37 


0.58 


0.063 


G1588 


4 


4 


0.73 


0.11 


0.94 


0.026 


0.88 


0.067 


0.63 


0.17 


0.93 


0.036 


0.73 


0.15 


G1667C 


8 


7 


-0.31 


0.23 


0.89 


0.0016 


-0.04 


0.40 


0.11 


0.40 


0.35 


0.23 


-0.15 


0.34 


G1674 


4 


5 


0.35 


0.36 


0.91 


0.045 


0.60 


0.23 


0.59 


0.18 


0.81 


0.036 


0.85 


0.028 


G1680 


4 


4 


0.96 


0.022 


0.85 


0.085 


0.89 


0.087 


0.95 


0.021 


0.13 


0.39 


-0.04 


0.50 


G1699 


4 


4 


0.57 


0.21 


0.43 


0.32 


0.25 


0.38 


0.32 


0.35 


0.89 


0.062 


0.48 


0.27 


G1832 


6 


5 


0.64 


0.083 


0.97 


0.0006 


0.05 


0.47 


-0.17 


0.37 


0.49 


0.15 


-0.02 


0.49 


G1849 


6 


6 


0.80 


0.063 


0.71 


0.069 


0.44 


0.24 


0.71 


0.077 


0.82 


0.041 


0.84 


0.025 


G1876 


4 


4 


0.76 


0.14 


1.00 


0.0005 


0.58 


0.24 


0.88 


0.054 


0.88 


0.053 


0.93 


0.032 


G1877 


<4 


5 














0.90 


0.032 


0.57 


0.20 


0.79 


0.088 


G1887 


4 


4 


0.85 


0.12 


0.91 


0.063 


0.63 


0.21 


0.98 


0.018 


0.82 


0.10 


0.86 


0.082 


G1908 


6 


8 


-0.24 


0.33 


0.89 


0.016 


0.31 


0.30 


0.12 


0.39 


-0.16 


0.37 


0.07 


0.44 


HIP12961 


6 


6 


0.35 


0.25 


0.30 


0.32 


0.63 


0.098 


-0.38 


0.23 


0.83 


0.018 


-0.78 


0.033 


HIP19394 


<4 


6 














0.39 


0.21 


0.08 


0.45 


0.26 


0.31 


HIP38594 


<4 


4 














0.59 


0.19 


0.19 


0.41 


-0.15 


0.43 


HIP85647 


5 


5 


0.94 


0.011 


0.99 


0.0016 


0.94 


0.012 


0.92 


0.019 


0.88 


0.037 


0.77 


0.087 



ity level (iMeunier & Delfossell2009h . This will produce the ob- 
served positive correlation between the two indices for higher 
values of S . We could also argue that this is an effect of the dif- 
ficulty of obtaining statistically significant correlations for stars 
with lower activity values but this seems not to be the case be- 
cause we found significant correlations for the cases of Na i and 
He I for low values of S (middle and right panels of Figure [S]). 
We therefore attribute this trend to the decrease of the impor- 
tance of filaments and the increase of the contribution of plages 
to the Hof index as the activity level of the stars increases. 

We found no dependence of the correlation between these 
indices and color. We note, however, that this was observed with 
very low signal-to-noise ratio for the (bluer) spectral regions that 
contain the Ca ii H and K lines. 

6.2. 5 Call vsNai 

Table m Figure [8] (middle panel), and the time-series presented 
in Figures [2l|7l show that the Nai behaves similarly to the 5'caii 
index^ 

iDiaz et al.l (l2007ah studied the Nai Dl and D2 lines for a 
sample of late-F to mid-M and found that these lines could be 
used as chromospheric activity indicators for very active stars. 
In our sample of low-activity stars, 16 out of 23 (around 70% of 
the sample) have a strong correlation (with FAP < 0.05) between 
the Na l index and our S -index. We have no negative values of 
the correlation coefficient and all the stars except three have p > 
0.5. Note that neither of these three stars passed the variability 
test for these two indices. Furthermore, we found no relation 



between the correlation coefficient and the average values of S , 
Na I, or with V - 1. This means that the Na I lines can be used 
as a proxy of the long-term activity as measured by Call for 
low-activity early-M dwarfs. This can be important because the 
Can H & K lines are located in the blue part of the spectrum 
and in these late-type stars the flux in this region is very low 
compared to the redder regions wh ere the Nai Dl and D2 lines 
are located. That lDiaz et al I (l2007ah have not found this trend for 
the least active stars of their sample might be because they used 
medium-resolution spectra, while we used spectra obtained with 
high resolution. 

6.3. ^caii v'sHei 

When we studied the variability of the activity indices in Section 
m we noted that the He l index was not showing as much vari- 
ability as the other three indices. Only three stars present a vari- 
ability with P- value below 0.05. Now, when we studied the cor- 
relation coefficient between this index and the 5'caii, we ob- 
served that the correlation between the two is very weak. Only 
three stars show correlation below a FAP value of 0.05 and the 
correlation coefficients vary between -0.50 and 0.94. This fur- 
ther confirms that the He I index is a poor activity proxy (as 
measured by the 5" -index) for early-M dwarfs. Figure [8] (right 
panel) shows the variation of the coefficient with mean 5" -index. 
Although there is a tendency for positive correlations as shown 
in the other indices, the coefficients are more scattered than for 
the 5'-Nal case and less significant. A similar trend as the one 
observed for Ha can also be observed here, where for stars with 
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Fig. 8. Pearson correlation coefficients for the relations between 5'caii, Ha, Nai, and He i as a function of mean 5'caii for the 23 stars 
with values of S index. Open squares are values with FAP < 0.01, open circles for 0.01 < FAP < 0.05, and "plus" ("+") symbols 
for FAP> 0.05. 



lower < 5'caii > there is a tendency for lower or negative values 
of the coefficient. 

6.4. Ua vs Na i, Ua vs He i, and Na i us He i 

Because the correlation between the calcium and the sodium 
lines is very strong for almost all stars in our sample, it is ex- 
pected that the relation between these two indices and Ha will 
be similar. This is the case, because the majority of the stars that 
show strong correlations here also have strong correlations be- 
tween S and the hydrogen line, except for G1581, which shows 
an anti-correlation for the two cases. Therefore, the discussion 
in Section [6J1 could also make sense in this case. 

Nine stars out of 30 show strong a correlation between Ha 
and He l (around 30% of the sample) with four more stars having 
0.05 < FAP < 0. 1 . This shows that the He l index is more similar 
in behavior to Ha than to 5'caii, where only three out of 23 stars 
showed a strong correlation. 

As expected from the S-Re l comparison, the correlation be- 
tween the Nal and He l indices is weak. Only 17% of the stars 
have FAP values lower than 0.05. 

6.5. Nightly averaged correlations 

Although the aim of this work was not to study the correlations 
between the four indices on short timespans (our concern was 
with long-term, cycle-type variations), we also calculated the 
correlation coefficients for the relations S caii-Ha and S caii-Na i 
for the nightly averaged data. The correlations observed for the 
binned data were maintained. The trend observed for the cor- 
relation between S and Ha when plotted against average S is 
also present and Na i maintains a very good correlation with the 
S index that is independent of the activity level. The FAP val- 
ues in general were lower because here we had more data points 
and therefore it was more difficult to obtain higher correlations 
by permutating individual points. Again, this result supports the 
use of the Na i index as a proxy of activity of M dwarfs even for 
high-frequency variations as does the fact that the Ha activity 
indicator does not vary linearly with activity as measured by the 
S index. 

This study will be presented in more detail in a future paper. 

6.6. Mean activity indices and color 

We compared the average values of the four activity indicators of 
each star to look for correlations (Figures [9] and [TOl). We found 
that the mean values of 5'caii, Nai and Hei are linearly corre- 
lated, but the same is not true when we compared them with the 



hydrogen line. Because these activity proxies were not corrected 
for the contribution of photospheric flux, we expect these indices 
to be dependent on the spectral type when we compare stars with 
difl'erent colors. As we can see in Figure [TT] the same three in- 
dices that have a linear correlation between their average values 
appear to have a similar trend with V - I: all of them decrease 
with color. On the other hand, the Ha index increases with V-I. 

We corrected the average values of the four indices using the 
residuals of the best-fit lines shown in Figure [TT] Figures [12] and 
[T3l show the relations between the residuals of the indices after 
correction. The linear trends between < Na I > vs < 5" > and 
< He I > vs < 5" > were maintained with similar correlation co- 
efficient values except for < He i > vs < Na i >, which decreased 
p from 0.80 to 0.55. 

Although the average values of these three indices are corre- 
lated, the scatter is still high. It is known that the fraction of ac- 
tive M dwarfs is not constant with spectral type but, increases for 
spectral types later than M4 (iDelfosse et al.lll998k . This could 
be a source of biasing toward the redder end of spectral color. 
However, in this sample only two stars have spectral type later 
than or equal to M4 and the redder of them (G155 1 , with spectral 
type M5.5) was discarded from these corrections. 

It is clear from Figures [12] and [13] that direct comparisons 
between Ha and the other three indices are not possible because 
there is no correlation between them. This reinforces the idea 
that the hydrogen line does not measure activity in the same 
way as Ca ll, Na l, or He l for early-M dwarfs. The strong intra- 
star correlations observed between the S and Na i lines are also 
observed when the activity level of diff'erent stars is compared. 
Therefore, the sodium lines can be used not only to measure the 
relative activity variations for a given star, but also to compare 
activity between stars. 



7. Conclusions 

We used more than five years of HARPS high-resolution spectra 
of a sample of M0-M5.5 stars from the HARPS M-dwarf planet 
search survey to study the behavior of the 5'caii, Ha, Nai, and 
He I chromospheric activity indicators. The data were binned to 
average-out unwanted short time- scale variations. We first ad- 
dressed the question of the mid- to long-term variability of the 
activity indices using F-tests. We found that some stars in our 
sample were showing obvious signs of these variations in 5'caii, 
Ha, and Na i. The same was not as evident for the He i index. 
The most obvious cases of detected long-term activity variations 
are GU, G1273, G1433, G1436, G1581, G1588, G1667C, G1832, 
G1849, G1877, G1908, and HIP85647. This selection was made 
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by choosing stars with statistically significant variation in at least 
two indices. 

We also performed F-tests to detect possible maxima or min- 
ima of activity. GU and G1433 have statistically significant max- 
ima in at least two indices. G1581 appears to be in a statisti- 
cally significant minimum of activity in all indices except for 
Ha, which shows a maximum. 

Our next step was to compare the activity indicators using 
th e Pearson correlation coeffi cient. Similarly to what was found 
by lCincunegui et all (l2007bh we obtained a great variety of cor- 
relations between the 5'caii and Ha in the range -0.8 < p < 1. 
We found that there is a trend of the correlation coefficient with 
the level of activity as measured by the S -index. For the least 
active stars there is a tendency for low- or anti-correlation, while 
for the most active stars in our sample the correlation becomes 
strongly positive. 



We attribute this effect to the fact that as activity increases, 
the Ha lines first start as a weak absorbing line, then become 
strongly absorbing, and as the activity further increases, they 
start to fill in, becoming weak absorbers again, and finally at 
hi gh-activity levels become a full emission line, as pr edicted 
by ICram & Muffanl (II l79b ( s ee also TCram & Giampapal (Il987h 
IStauff'er & HartmanJ^l), IWalkowicz & HawlevI (l2Q09l) ). A 
relation like this would produce the correlation observed in 
Figure [8] (left panel). This ambiguous relation between the cal- 
cium and Ha lines was also observed in Figure [T] where we 
compared the shape of the activity lines for different activity 
levels. This behavior might be a consequence of the different 
contribution of filaments to both indices, their different spa- 
tial distribution in the stellar disk, and the fact that its con- 
tributio n to the Ha index saturate s at a given level of activ- 
ity (see iMeunier & Pelf ossd [2009b . We thus conclude that, al- 
though the Ha is able to detect long-term activity variations in 
early-M dwarfs, it is not the same physical phenomena that is 
being probed by both indices and thus their variations cannot 
be naively compared. Nevertheless, the simultaneous study of 
^Caii (or Nai) and Ha might be useful to gather new informa- 
tion about the presence of different activity-related features such 
as filaments and plages in stellar disks. 

The He i index has not only the smallest variation of the three 
lines, we also found it to correlate less well with activity and thus 
we do not recommend its use for the long-term activity study of 
M dwarfs. 

We found a very good correlation between 5'caii and Nai. A 
great number of stars display a similar behavior for both indices 
and there is no relation between the correlation coefficient and 
the activity level. This index consequently is a proxy of activity 
variations as measured by 5'caii- Because these stars have more 
flux in the redder part of the spectrum, the Nai index should 
be preferred over the Ca ii lines because the measurements will 
have a far better signal-to-noise ratio. We can then extend the 



16 



J. Gomes da Silva et aL: Long-term magnetic activity of a sample of M-dwarf stars from the HARPS program 



p = 0'.33 ' ' ' ' 


' ^ = 23 


- ' ffl ^ 


" 1*1 1*1 _ 
* - 



O-C < Scall > 



0.05 





0.08 




0.06 




0.04 


A 










0.02 


V 




o 


0.00 


o 






-0.02 




-0.04 




-0.06 



p = 0.77 






= 23 
















^ " " 




















* 















D.0020 
D.0015 
3.0010 
D.0005 
D.OOOO 
3.0005 
D.OOlO 
0.0015 



p = 0'.71 ' ' ' ' 


^ = 23 


^^^^ 








1 a 1 





-0.04 -0.03 -0.02 -0.01 0.00 0.01 0.02 0.03 0.04 0.05 
O-C < Scall > 



-0.04 -0.03 -0.02 -0.01 0.00 0.01 
O-C < Scall 



0.02 0.03 0.04 0.05 



Fig. 12. Comparison between the average values of 5'caib Ha , Na i, and He i for the 23 stars with values of S index after correction 
for color. Symbols as in Figure [9l 




-0.02 0.00 0.02 0.04 
O-C <NaI> 



0.012 
0.010 
0.008 
0.006 
0.004 
0.002 
0.000 
-0.002 
-0.004 
-0.006 



p = O.'lS 




1 ' i 1 ' 


^ = 29 




* 



































-0.0020-0.0015-0.0010-0.00050.0000 0.0005 0.0010 0.0015 0.0020 
O-C <HeI> 





0.0020 




0.0015 




0.0010 


A 




Oi 


0.0005 


K 




V 


0.0000 


o 






-0.0005 


o 






-0.0010 




-0.0015 




-0.0020 





^ = 29 













0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.0 
O-C <NaI> 



Fig. 13. Comparison between the average values of Ha , Na i, and He i for 29 stars after correction for color. G1551 was excluded 
from these plots because of its very high activity values. Symbols as in Figure [9j 



iDiaz et aL I (l2007ah suggestion of the use of these lines to mea- 
sure mid- to long-term activity in the least active stars of later 
spectral type. 

In a following paper we will compare the behavior of the 
long-term activity from these results with the radial-velocity 
measurements and the cross-correlation function parameters bi- 
sector span, full-width-at-half-maximum, and contrast, with the 
aim of investigating if long-term cycle-like activity variations are 
able to influence the measured radial velocity of M dwarfs. 
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